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Abstract 

We introduce in this paper a new tool to prove the convergence 
of the Overlapping Optimized Schwarz Methods with multisubdo- 
mains. The technique is based on some estimates of the errors on 
the boundaries of the overlapping strips. Our guiding example is an 
n-Dimensional Linear Parabolic Equation. 

1 Introduction 

In the pioneer work [9], ^0\, [H], P. L. Lions laid the foundations of 
the modern theory of Schwarz Methods. With the development of parallel 
computers, the interest in Schwarz Methods have grown rapidly, as these 
methods lead to inherently parallel algorithms. However, with Classical 
Schwarz Methods, high frequency components converge very fast, while low 
frequency components converge slowly and that slows down the performance 
of the methods. By replacing Dirichlet Transmission Condition in Classical 
Schwarz Methods by Robin or higher order Transmission Conditions, we can 
correct this weekness of Schwarz Method. The new methods are called Opti- 
mized Schwarz Methods and have been introduced in [5], [6]. Since then, the 
convergence properties of the Optimized Schwarz Methods have been stud- 
ied deeply, based on the following two main tools: Energy Estimates and 



1 



Laplace and Fourier Transforms. Energy Estimates allow us to study the 
convergence of the methods in the case of nonoverlapping subdomains. With 
Energy Estimates, both linear and nonlinear problems have been studied and 
Optimized Schwarz Methods have been proven to converge, while applying 
to these equations (see for example, the papers [1], [7], |8]). On the other 
hand, Laplace and Fourier Transforms allow us to study the convergence 
of the Overlapping Optimized Schwarz Methods, but for only a few simple 
equations (see, for example [2], [5], [6]), and the convergence problem 
of the Domain Decomposition Methods with Robin Transmission Conditions 
still remains an open problem up to now. 

In this paper, we introduce a new tool to prove the convergence of the 
Optimized Schwarz Methods for multisubdomains and apply it into an n- 
dimensional linear parabolic equation of the following form 

The idea of the technique is to estimate carefully the difference between the 
values of the errors at the boundaries of the overlapping strips. The technique 
has the potential to be applied to many other kinds of Partial Differential 
Equations including nonlinear ones. 



2 Problem Description and Main Results 

We consider the following parabolic equation 

\ u{x,t) = g{x,t), on dn X {0,T), 
u{x,0) = g{x,0), on fi, 

(2.1) 

where Q = D x {a, (3), D is a. bounded and smooth enough domain in R"" ^. 
We impose the following conditions on the coefficients of (12.11) 
(Al) For all i,j in {1, . . . , /}, aij(t) = ajs{t). There exists > Q such that 
A{t) = (ajj(t)) > i^qI for all t belongs to (0,T) in the sense of symmetric 
positive definite matrices. 

(A2) The functions bounded in C°°(]R); / and g are bounded 
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functions in C^{n x (0,r)). 

With the conditions (Al) and (^2), Equation f l2.ip has a unique bounded 
solution u in C°°{{0,T) x Q). The proof of this result can be infered from 
Theorems 9 and 10, page 71 [3J. 

We now divide the domain Q into / subdomains, with Qi = D x {ai,bi) 
and a = ai < a2 < bi < ■ ■ ■ < ai < < bj = 13. The Optimized Schwarz 
Waveform Relaxation Algorithm solves / equations in I subdomains instead 
of solving directly the main problem (12. ip . The iterate j^k in the l-th domain, 
denoted by uf, is defined by 



^ = '-^^^+put^\;ai,-) in D x (0,T), 

+P^^f(-,^/,-) = '"'it'''" +K+i^(->^''-) inDx(0,T), 



k 



(2.2) 

here, p is a constant and for each vector x in M", we denote a; = (X, x„), 
with X G and x„ G M. Each iterate inherits the boundary conditions 

and the initial values of u: 

u'l{x,t) = g{x,t) on {dQj n dQ) x (0,T), nf(a;,0) = g{x,0) in 

and a special treatment for the extreme subdomains. 



A bounded initial guess in C°°(f2 x (0,T)) is provided, i.e. we solve at 
step Equations (12. 2p . with boundary data on left and right 



dXn 

dul{-,bi,- 



+ pul{-,k,-) = h\;bi,-) mDx{0,T). 



By using an induction argument and the same arguments as in Theorem 2, 
page 144 [3], we can see that each subproblem (12. 2p in each iteration has a 
unique solution. Theorem 10, page 71 [3] shows that these solutions belong 
to C°°(f2 X (0,T)). This means that the algorithm is well-posed. 
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Denote by ef the difference between and u, and substract Equation 
with the main equation (I2.ip . we get the following equation on ef 



f ¥-E:,=i«MWSt + Sr=i^^WS + c(t)ef = infi,x(0,r), 
+ Pefi; au-) = '-^^ + pet,\;a, ■) in D x (0, T), 

^ +Pef = ^^^^ in D x (0,r). 

(2.3) 

Similarly, each iterate inherits the boundary conditions and the initial values 
of u 

ef = on {dni n on) x (O, r), ef (-, 0) = in ni, 
and the special treatment for the extreme subdomains, 

«t(-,a,-) = 0, u%,/3,-) = 0. 
The following theorem states that the algorithm converges. 
Theorem 2.1 Let be a strictly positive function in C^(R) such that 

— max — [Xr, 

I-nGM \ if 

is large enough, the Optimized Schwarz Waveform Relaxation Method con- 
verges in the following sense 



lim max 

A;^oo/G{l,...,/} 



— |- exp(px„)j (/j(t) 



0. 



C(n,x(0,T)) 



Moreover, for I in {1, . . . , /}, the sequence {wf — u} converges pointwisely to 
as k tends to infinity. 

Remark 2.1 We can see that if we choose (f{xn) = exp{—^Xn), then if 'y 
is large enough, — max^,^g]R ^^(a^n)^ is large enough. The condition of our 
theorem is then satisfied. 

Remark 2.2 Since aij, hi are functions of t, and the domain is divided in 
to n-subdomains, we cannot use Fourier and Laplace Transforms. Moreover, 
since the subdomains are overlapping, the Energy Estimates Method cannot 
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be used in our case. In the next section, we introduce a new technique to prove 
the convergence of the algorithm, the technique is based on the observation 
that we can estimate the difference between the values of ef on the boundary 
and in the interior. 

Remark 2.3 The result in the theorem remains true if we let aij, hi be 
bounded and continuous functions of t and x, but not depend on the n-th 
space variable x„, as we can see in the proof in the following section. 



3 The Convergence of the Algorithm 

This section is devoted to the proof of Theorem 12.11 We divide the proof 
into two steps. 

Step 1: The Error Estimates. 

For /c G N and i G {1,...,/}, setting ef = efexp(px„), we get ef = 
ef exp(— px^). Equation (12. 3 p then leads to 



+ (c - pbn - p^an,n)ef = 0, in x (0, T), 
^^^^mDx(0,n 



e 



on (dnjodn) X (o,r), 



L ef(-,-,0) = Oinfi 



'J' 



and for the extreme subdomains, 

e^(-,ai,-) = 0, e%,bj,-) = 0. 

Setting 



(3.1) 
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we infer from Equation f l3.ip that for / in {2, . . . , J — 1} 

+ (c - - p^anM = 0, in Qi x (0, T), 
i.f(-,a,,-) = /vf:-'(-,a;,-) inDx(0,T), 
z/f(-,6;,-) = <-i'(-,&i,-) inDx(0,T), 



on ((9fij nSfi) X (o,r), 



z/f(.,.,0) = in fij. 

(3.2) 

We can observe that these systems are with Dicrichlet tramission conditions. 



On X (0,T), we define 
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k \ 2 



where is a strictly positive function in C^(]R) to be chosen later, with the 
notice that — maXa-^gK {^{xrSj is large enough. Our purpose is to construct 
an operator £ of $, such that is negative and then on £, we can apply 
the maximum principle to get some estimates on the boundaries for $. With 

these estimates, we can direclty infer some good estimates for ^j- and that 
lead to our result on the convergence. We now consider the following operator 

(3.3) 

A simple calculation gives 

m) - _ ± _ ± (3.4) 

\ %,]=! * ■> / i,j=l * J 

Since the second term on the right hand side of the previous inequality is 
negative, it directly leads to 
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k\2 



■ ■ 1 *"^7 

i,j = l ■> 



- 5^ 2a,„0Vi^f ^ + - an,n^J M^if- (3.5) 

Our purpose is to transfer the right hand side of f l3.5p into the sum of a 
negative term and a term of in order to do that, we replace (13. 2 p into 
(13. 5 p and get the following bound for 

1=1 * " 

" r)u^ fin' th"\ 

- Xj2a,.0V^z'^ + (^^ - an,n^J (3.6) 
Replacing 

— = 20(/?i/f-^, forie 

OXi OXi 



and 



OXfi dxyi 

into (13. 6p . we get 

- 2. 5^ + ^2pan,. + fc. + a„,„- J ^ (3.7) 

i,j=l J \ ' / 



+ 2_, \ P(^i,n + Oi + ain — J 



dxi 



< <t>f{iyif ( — - ann— + fln.n-T - 2(c - p6„ - p^a,i,n) + (2pan,n + bn)^ 
\(p (f) (j)^ (f) 

We now get the formula for 2, 

£(*) := 5^ a.,— + (2pa,„ + 6„ + «„..-)— (3.8) 

l,j = l 

n-1 



then if we choose Lp such that — maXi.^giR {^{xnYj is large enough, since an^m 

bn, c, ^, ^ are all bounded in C(M), we can obtain a negative sign on the 
right hand side of fl3.8p . which means is negative. 

Since < 0, the maximum of $ can only be attained on the boundary 
of X (0, T). Since $ = on dQiCid^ and on i7 x {0}, we have the following 
three estimates. 

Estimate 1: I < I < I. 

The maximum value(s) of $ can be achived on both D x {a/} x [0, T] and 
D X {bi} X [0,T] and 

(z.f(X,x„,t))20(x„)(/.(t) < (3.9) 
< max{ max {(i/f (X, a,, t))='0(a;)(/?(t)}, max {(i^f (X, 6;, 

Dx[0,T] Dx[0,T] 



Estimate 2: 1 = 1. 

The maximum value(s) of $ can be achived on both D x {ai} x [0,T] 
and D x {61} x [0, T]. If the maximum of $ is achived on D x {ai} x [0, T], 

d± 
dn 

compute 



then at the maximum point, we need that > due to Hopf 's Lemma. We 



Since 



we get 



dn 



ai,t) 



dx 

dxl dx„ 



4>{ai)(p{t) - 



(P'ia.Mt) 



-0(ai)(/?(t) 



dxldxn \dxnj 0(ai) 



'dt 



n 



— h > (pai,„ + bi)- h (2pa„,„ + &„)- — 

An,n)et = 0, inr^i X (0,T), 



1=1 
+ (c - p6„ 



n-1 



ai, 



+ {2pan,n + bn)-g^{.,ai,.) + {c-pbn-p^an,n)e^{.,ai,.) = 0, inDx (0,T). 

Combining this equation and the fact that ef(.,ai, .) — on D x {0,T), we 
can deduce 



dxl 



and as a consequence, we can write |^ in a different way 



dn 



.,ai, .J 



-(l){ai)^{t) 



del 
dxr, 



2j9 + 



£(«i) 

0(ai) 



With the functions satisfying 



2p + 



+ 



we can see that 



>o, 



dn 



.,ai,.) < 0; 



which means that the maximum of $ can be achived only on D x {hi} x [0, T] , 
then 

{v\{X,Xn,t)f4>Mv{t) < max {{iy',{X,bi,t))'(f>{h)ipit)}. (3.10) 

Llx[0,T] 



Estimate 3: I ^ I. 

The maximum value(s) of $ can be achived on both D x {a/} x [0, T] and 
D X {67} X [0, T]. If the maximum of $ is achived on D x {67} x [0, T], then 
at the maximum point, we need that > due to Hopf's Lemma. Similar 
as in Estimate 2, we can get 



dn 



.,bi,t) 



dXn 

= <P{bjMt) 



+ 



d^el del , [ deiy (j)'{h_ 
dx"^ dxn 

de\ ^ ^ 



2p + 



<^(^>7) 
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With the functions satisfying 



bJt) \ (p'ihi) 



an,n{t)J (f){bi) 
we can see that 

which means the maximum of $ can be achived only on D x {a/} x [0,r], 
then 

iuXX,Xn,t)f(p{xnMt)< max {{iyXX,ai,t)f<P{ai)ip{t)}. (3.11) 

Dx[0,T] 



Step 2: Proof of convergence, 



hm max 

fc->oo «g{1, ...,/} 



0. 



c{nix{o,T)) 



In the proof of convergence, we will use the three estimates fl3.9p . (13.101) 
and (13.111) by fixing the function {p{t) and replacing </)(x„) by appropriate 
functions (j)^, (pi, (j)^, (p^ {i E {1, . . . ,/}) in each subdomain. Before coming 
to the details of the proof, we need a notation 

= ^^max ^ ||(z/f)20(^||^(^^^^^^. 

Step 2.1: Estimate of the right boudaries of the sub-domains. 
Consider the I-th domain, at the k-th step, (13.111) infers 

{iy'}{X,Xn,t)y'^j{xnMt) < max {{iy^{X,ai,t)yMai)ip{t)}, 

Dx[0,T] 

where is a strictly postive function and will be chosen later. 
Replace Xn by we get 

{uXX,bj.i,t)fMbj.iMt) < max {{i^XX,aj,t))^j{ajMt)}. 

Dx[0,T] 

Since ly^iX, bj^,, t) = ly'l+HX, bj^,, t), 

{u^j^l{XM-iM'4>i{h-Mt) < max {{u'}{X, aj,t)yMaiMt)}. 

Dx[0,T] 
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The inequality becomes 

iu^^l{X,bj^,,t)fcp{t) < max a,, t)) VW}- 

We can choose 67 such that -^'.[""'K < 1, and deduce 

tt>i{bl-i) 

iu^+l{X,bj-i,t))Mt) < T7r\^^- ^3.12) 
Moreover, on the (/ — l)-th domain, at the {k + l)-th step, (13.91) leads to 

max{ max {{u'}+l{X,bi^i,t))%_^{bi-i)ip{t)}, 

Dx[0,T] 

max {{u^+l{X,ai.i,t))%_,{aj^,)^{t)}}, 

Dx[0,T] 

where 4'j_i is a strictly postive function and will be chosen later. 
Since z/,^+i(X,6,_2,t) = u^+^iX,bj.2,t), 

{iy^^',{X,bj^2,t)m^,{bi.,Mt)< 

max{ max {{u'}+l{X,bi^i,t))^'^j_,{bj^i)ip{t)}, 

Dx[0,T] 

max {{iy'}^l{X,ai.,,t)f^j_,{ar_,Mt)}}. 

Dx[0,T] 



Hence 



iu^liix,bj^,,t)yvit)< 

maxima max {(1.^+^ (X, t)) }, 

max ||-4r4(.t_+^(X,a,_,,t))VW 
Combining this inequality with fl3.12p . we get 



{ur^iX, 6,_2, t))Mt) < max { tl^l^M^E,, tl^li^E, 
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fc+1 



We choose such that 

^j-i(^j-i) 0/(aj) ^ 0j-i(aj-i) ^ ^ 

^/_l(f'/-2) 0/_i(b/-2) 

which is equivalent to 

As a consequence, 

(z.t+'(^,&/-2,t))VW < m^max{i?,,i?,+i}. (3.13) 

Using the same techniques as the ones we use to achive f l3.12p and ( I3.13P , 
we can prove that 

^^^■(X,6,_„t))VW < %^ii^^^=^max{E,,...,i?,+,_i}, (3.14) 
where 4>i-j+i is a strictly positive function satisfying 



<1, 



withi = {l,...,/-l}. 

Now, with (13.1 op . we can choose a strictly postive function 0^ such that 
0*(^/) > 0*(«/), then 

iu^iX,bj,t)f^M^it) < max {{u'^iX,ai,t)f^MMt)}, 

Dx[0,T] 

and clS cl result 

{iy^{X,br,t)YMbiMt) < max {(i.ti'(X, a,, t))20,(a,)(/p(t)}, 

Dx[0,T] 



which implies 



[uHX,bj,t))Mt)<^^E,_,. (3.15) 
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Step 2.2: Estimate of the left boundaries of the sub- domains. 
Consider the 1-th domain, at the k-th step, fl3.1ip infers 

{iy1{X,Xr„t)fMxnMt)< max {{u'l{X,b,,t)fMhMt)}, 

Dx[0,T] 

where 01 is a strictly postive function and will be chosen later. 
Replacing Xn by 02, we get 

K^(X,a2,t))VW < -^^^ {(z.f(X,6i,t))V(t)}. 

(pl{a2) Dx[0,T] 

Since z/f(X, a2,t) = h'2~^^{X,a2,t), then 

(z.^i(X,a2,t))VW < max {(j/f (X, 61, t)) VW}- 

(/>l(a2) Dx[0,T] 



We choose such that 

MK 



< 1, 



01(02) 

and deduce 

(z.^i(X,a2,t))VW < TT^^;^- (3.16) 
Moreover, on the 2-th domain, at the {k + l)-th step, fl3.1ip leads to 

{ul^+\x,x^,tyfMxMt)< 

max{ max {{u^+\X,b2,t))^Mb2Mt)}, 

Dx[0,T] 

max {{iyt\X,a2,t)fUa2Mt)}}, 

Dx[0,T] 

where 02 is a strictly postive function and will be chosen later. 
Since z/2"^^(X, 03, t) = z/3+^(X, 03, t), then 

{u',+\X,a,,t))%{a,Mt)< 

max{ max {{u^+\X,b2,t)YMb2Mt)}, 

Dx[0,T] 

max {{u^+\X,a2,t))'Ma2Mt)}}- 

Dx[0,T] 
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Hence 

iu',+\x,as,t)rm< 

max 114^ max {{ul+\X,h,t))'^{t)} 

[Mas) D^T] 

ja^x{iu',^\X,a,,t))Mt)} 

02 (as) Dx[0,T] 

Combining this with (]3.16p . we get 

(u^^X, a3, t))Mt) < max | t^E,^,, Ma^kib^E, ^ . 

[Mas) 02(03) 01 (^2) 

We choose 02 such that 

02(^2) _ 02(02) 0i(^i) ^ ^ 
02(03) 02(03)01(02) 

which is equivalent to 

02(^2) _ 0l(fel) ^ ^ 
02(02) 01(02) 

We then obtain 

{lyt^X, as, t))Mt) < max{i?,, E,^,}. (3.17) 

02(03) 

Using the same techniques as the ones that we use to achive (13.161) and 
(I3.17p . we can prove that 

(/.f ^■-i(X,a„t))V(t) < \'^^rl^ me.x{E,, . . .,E,+,.,}, (3.18) 

0i-i(Oj) 

where (pj-i is a stricly positive function satisfying 

4>j-i{bj-i) ^ ^ 

^j-l(Oj) 
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with j belongs to {1, . . . , J — 1}. 

Now, with fl3.11l) . we can choose a strictly positive function 0^, such that 
<j)*{bi) < ^*(ai), and get 

{ly',{X,a^,t)fMalMt) < max {{u',{X,h,t)fMbiMt)}, 

Dx[0,T] 



which is equivalent to 



(z/f(X,ai,t))2(/),(ai)(^(t) < max {{ut\X,b,,t)yMbiMt)}. 

Dx[0,T] 



That implies 



[u'}iX,a,,t))Mt)<i^E,_,. 



(3.19) 



Step 2.3: Convergence result. 

From daUD, fl37[5|l . f l318|) and fl319|l . there exists 7 in (0, 1) such that 
(z/f+^(X,a,,t))V(t) <7niax{Efc,...,£:fc+,_i}, for / G {1, . . . , /}, (3.20) 



and 



(//f+^(X, t))V(t) < 7max{^fe, . . . , ^^+7-1}, for / G {1, . . . , /}. (3.21) 
Using (13. 9p for = 1, we have that 

(i^f(X,x„,t))V(t) < (3.22) 
< max{ max {(//f (X, a^, t)) V(t)}, max {(z/f (X, 6^, t)) V(t)}}- 

_Dx[0,r] Dx[0,T] 



Combining (Km . (Km and ([322]), we get 

Ek+i < 7max{Efc, . . . 
Hence, Ek tends to as A; tends to infinity. That gives 



(3.23) 



lim max 
k^oo ie{i, ...,/} 



del 

dXr 



k\ 2 



0. 



c{nix(o,T)) 

Step 3: Proof of convergence: for / in {1,...,/}, the sequence {ef} 
converges pointwisely to as A; tends to infinity. 
Since for / in {1, ... , J}, 



lim 



<9ef 



dXr. 



c{nixio,T)) 
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then I i^^j '^{t) 1 converges to pointwisely and the sequence is bounded 

by a constant Mq. Since (p is strictly positive on [0,T], there exist posi- 
tive constants Mi, M2 such that Mi < ip < M2. That means the sequence 

I dx~ I converges to pointwisely and is bounded by a constant M3. 

For I — 1, with a fixed value of {X,t), the Lebesgue Dominated Conver- 
gence Theorem confirms that for Xn belongs to [ai, bi], 

converges to as tends to infinity Hence the sequence 

{e'iiX,Xn,t)-e1iX,ai,t)} 

converges to as /c tends to infinity Since e^{X, ai,t) = 0, then the sequence 
{e^} converges to pointwisely 

For I — 2, with a fixed value of {X, t), again by the Lebesgue Dominated 
Convergence Theorem, for Xn belongs to [02,62], the sequence 

converges to as tends to infinity Hence the sequence 

{e^(X,x„,t) -e^(X,a2,t)} 
converges to as A; tends to infinity Since 

^{X, 02, t) + pe'^iX, a2, t) = -^(X, 02, t) + pet\X, 02, t), 

and the sequences {e^}, < \q^\ \ converge to pointwisely for / belongs 

to {1, . . .,/}, we can deduce that e2{X,Xn,t) converges to as A; tends to 
infinity. 

By similar processes, we can prove that for / in {1, . . . , /}, the sequence 
{ef} converges pointwisely to as A; tends to infinity. This concludes the 
proof. 
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